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Time-dependent studies of inclusive charm decays to multibody self-conjugate final states can be 
used to determine the indirect CP-violating observable Ar and the mixing observable yep, provided 
that the fractional CP-even content of the final state, P+, is known. This approach can yield 
significantly improved sensitivity compared with the conventional method that relies on decays to 
CP eigenstates. In particular, D —» 7r + 7r _ 7r° appears to be an especially powerful channel, given 
its relatively large branching fraction and the high value of P+ that has recently been measured at 
charm threshold. 
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It is of great interest to search for effects of indirect 
CP violation in time-dependent studies of neutral charm- 
meson decays. In the Standard Model indirect CP vi¬ 
olation is expected to be well below the current level 
of experimental precision (1), but many models of New 
Physics predict enhancements [ 2 ). A very important CP- 
violating observable is Ar, which is measured from the 
difference in lifetimes of the decays of D° and D° mesons 
to a CP eigenstate. In this paper it is shown how in¬ 
clusive self-conjugate multibody decays that are not CP 
eigenstates can also be harnessed for the measurement of 
Ap, provided that their fractional CP-even content, P+, 
is known. This new approach has the potential to im¬ 
prove significantly the knowledge of Ap and has become 
possible thanks to measurements of P+ that have recently 
begun to emerge from analyses of coherent charm-meson 
pairs produced at the ^(3770) resonance @,[i]. Further¬ 
more, it is explained how exploiting these decays can 
also provide a corresponding improvement in the preci¬ 
sion on yep , which is an important observable that de¬ 
scribes D°D° oscillations. For the purpose of concrete¬ 
ness the discussion is presented for the example decay 
D —» 7r + 7r - 7T°, although the results are valid for all self¬ 
conjugate multibody modes. Here and throughout the 
discussion D indicates a neutral charm meson; this no¬ 
tation is used when it is either unnecessary or not mean¬ 
ingful to specify a flavour eigenstate. 

Measurements with CP eigenstates 

In the D-meson system the mass eigenstates, £> 12 , are 
related to the flavour eigenstates D° and D° as follows: 

\D h2 )=p\D°)±q\D°), (1) 

where the coefficients satisfy \p\ 2 + \q\ 2 = 1 and 

r C pe^ CP = 1 (2) 

P 

The phase convention CP|D°) = | D°) is adopted. Indi¬ 
rect CP violation occurs if rep ^ 1 and/or cfcp ^ 0- 


Charm mixing is conventionally parameterised by the 
quantities x and y, defined as 


Af\ — I/o Ti — r 2 

f ’ V = 2T 


(3) 


where Af 12 and r 12 are the mass and width of the two 
neutral meson mass eigenstates, and T the mean de¬ 
cay width of the mass eigenstates. In the chosen con¬ 
vention D\ is almost CP even. The average of cur¬ 
rently available measurements gives x = ( 0 . 41 ^ 0 + 5 )% 
and y = (0.63to;°g)% 0. 

Consider an environment where charm mesons are pro¬ 
duced incoherently, such as the LHC or an e + e _ 13- 
factory, and are observed through their decay into a CP 
eigenstate of eigenvalue rj cp- Time-dependent measure¬ 
ments allow the decay widths T and F to be determined 
for mesons produced in the D° and D° flavour states, 
respectively. From these quantities the CP-violating ob¬ 
servable Ar and mixing observable yep may be con¬ 
structed: 


, f-f 

Ar = rj cp- - 

r + r 


yep =vcp 



(4) 


Assuming x , y, (rep — l/rcp) and <f>cp to be small, 
and assuming direct CP violation to be negligible, it can 
be shown [b] that these observables have the following 
dependence on the underlying physics parameters 


Ap 


Vcp 
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Expressions that also allow for the contribution of di¬ 
rect CP violation can be found in Ref. |7|. Thus in 
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the limit of CP conservation Hr vanishes and yep —> y. 
The average of currently available measurements, dom¬ 
inated by studies based on the CP-even eigenstates 
K + K~ and 7T+7T - , yields H r = (—0.058 ± 0.040)% and 
Vcp = (0.866 ± 0.155)% (Here the Hr average in¬ 
cludes new measurements from the LHCb Q and CDF [9] 
collaborations, in addition to the older set of results from 
LHCb [10J, BaBar 11] and Belle 12| that are considered 
in Ref. [if.) 

Introducing self-conjugate multibody decays and 
the CP-even fraction Fq. 

The CP content of an inclusive self-conjugate multi¬ 
body decay, for example D —> 7 r + 7 r _ 7 r°, can be measured 
with a sample of coherently produced DD pairs at the 
^(3770) resonance, such as that collected by the CLEO- 
c and BESIII experiments. A double-tag technique is 
employed in which one D meson is reconstructed in the 
signal decay of interest, and the other in its decay to 
a CP eigenstate. In such an event, and neglecting any 
CP violation, the quantum numbers of the ^(3770) me¬ 
son means that the CP eigenvalue of the signal decay 
is fixed. The CP-even fraction of the signal decay is 
given by F + = N + /(N + + N~), where N + (N~) des¬ 
ignates the number of decays tagged as CP-even (-odd), 
after correction for detector inefficiencies and the spe¬ 
cific branching fractions of the CP eigenstate tags em¬ 
ployed. In this manner F + has been measured for the 
decay D — > 7 r + 7 r _ 7 r° and found to be 0.973 ± 0.017, indi¬ 
cating the mode to be almost fully CP even [§, 4 }. 

Although CP violation is neglected in the currently 
available measurements of F + this assumption introduces 
negligible bias in the result. Both the Standard Model 
and theories of New Physics expect direct CP violation 
in charm decays to be < 10 -3 |l3|, a prediction which is 
compatible with existing experimental results [14]. Any 
effects will therefore be small alongside the measurement 
precision attainable with the CLEO-c and current BE¬ 
SIII data sets. Furthermore, the double-tag analyses 
performed at these experiments have no sensitivity to 
indirect CP violation at leading order in (x^y), as the 
DD system is produced at rest. For the specific case of 
D —> 7 T + 7 r _ 7 r°, a recent time-integrated high precision 
analysis by LHCb has revealed no evidence of any direct 
CP- violating effects [HI]. 

There is a simple relationship between F + and the 
parameters that describe the intensity and strong-phase 
variation over the phase space of the decay. The am¬ 
plitude of a multibody decay such as D —> is 

dependent on the final-state kinematics, which can be 
uniquely defined by the Dalitz plot coordinates S 12 = 
? 7 i 2 ( 7 r + 7 T°) and S 13 = m 2 ( 7 r _ 7 r°). The amplitude of a D° 
decay to a specific final state is given by Hdo(si2, S13) = 
ai 2 ,i 3 e*' 5l2 ’ 13 ) where the integral of |Hdo(si 2 , S 13)| 2 over 
the full Dalitz plot is normalised to unity. Consider the 
situation where the Dalitz plot is divided into two bins by 


the line S 12 = S 13 . The bin for which S 12 > S 13 is labelled 
— 1 and the opposite bin is labelled +1. The parameter 
Kj (A'j) is the flavour-tagged fractional intensity, being 
the proportion of decays to fall in bin i in the case that 
the mother particle is known to be a D° (D°) meson: 


Ki = / |ai2,13| 2 dsi 2 dsi 3 . 


(7) 


The parameter Cj is the cosine of the strong-phase differ¬ 
ence between D° and D° decays averaged in bin i and 
weighted by the absolute decay rate: 


l 


Ql2,13^12,13 COs(£i2,l3 — £ 12 , 13 ) 

\Zk,k, 


dsi2 dsi 3 . (8) 


A parameter Si is defined in an analogous manner for the 
sine of the strong-phase difference. 

The CP-tagged populations of these bins, N±, nor¬ 
malised by the corresponding single CP-tag yields, is 
given by jl 6 | 


N± = h D (Ki ± 2 a^KiKi + K^. (9) 

Here hp> is a normalisation factor independent of bin 
number and CP tag. When there is no direct CP 
violation in the decay Hpo(si 2 , S13) = ai 2 ,i 3 e *' 512,13 = 
ai 3 ,i 2 e *' 513,12 and so 

Ki = K-i, a = c-i and s t = — s-%. ( 10 ) 

Under this assumption, and the identities N ± = ]TL N±, 
and ]Tb K. ir = 1, it follows that in the two-bin case 

F+ = i (l + 2c lV / iCAU7) • (11) 


Measurements with inclusive self-conjugate 
multibody decays 

Now consider, for an incoherently produced D meson, 
the time dependence of a self-conjugate multibody decay. 
The time evolution of the Z?° to the point (si 2 ,si 3 ) is 
given by 

Apo{t, S12, S13) = ai2,i3 e lSl 2 , 13 g+(t) 

+ rcpe i<t,cp a^ 12 e iS ^g-{t), (12) 

where g±(t) = i[ e -*(M 1 - i r 1 / 2 ) t ± e -i(M 2 -ir 2 / 2 )tj_ Ignor _ 
ing terms of 0(x 2 ,y 2 ,xy) or higher, the rate of decay to 
that point is proportional to 

|Hi)o(t,Si2,si3)| 2 = e rt |ai2,i3 — 012,13^13,i2’"cprt x 
ycos(5i2,i3 ~ £13,12 — <!> cp ) + 
x sin(<5i2,13 — £13,12 — <t>cp) |- 

( 13 ) 
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Integrating this over the two bins of the full Dalitz plot With the inclusion of direct CP violation the expres- 
leads to the time-dependent decay probability sion for Ap fi becomes 


= [ \A D o(t, si2,si3)| 2 dsi 2 dsi 3 + 

J +1 

J |Ad° (t, S 12 , Sl3)| 2 dsi 2 dsi 3 (14) 
= exp(-Tt) 1 - r C p(2F+ - 1) x 


(y cos 4>cp ~ x sin <pc P )Tt 


(15) 


where use is made of the definitions of Ci , Si and the re¬ 
lations given in Eqs. m and CEB- Hence the width of 
the decay is approximated at first order in x and y by 


l + r C p(2F + -l) x (y cos (/) C p-x sin (j) C p) ■ (16) 


The time evolution for the D° decay to the point 
(si 2 ,si 3 ) is given by 


S12, S13) 


— e-^ OP ai2,i 3 e i5l2 - 13 fl-(t) 

TCP 

a 13A2 e i5 ^g+(t), (17) 


and thus the width of the D° decay is approximated by 


r « r 


1 + —(2F+ —1) x {ycosfycp+xsinfycp) . (18) 

TCP J 


Defining Ap ff = ¥- 4 ? and y^ P = — 1^ it follows 

that 

Ap ff « i( 2 F+ - \)y cos cf> C p ( r C p -— ) - 

2 V tcp J 

t:(2F+ - l)xsm.(j)cp (rep + ——) , (19) 
2 V r CPJ 

V e cp ~ t:(2F+ - l)y cos (j) C p (tcp + —) - 

2 V r GP J 

\{2F + - l)x sin (j) C p (r C p -— ) ■ (20) 

2 V r CPJ 

These expressions contain an additional dilution factor 
of (2F + — 1 ) in comparison to the CP-eigenstate relations 
of Eqs. d5j and © and are identical in the case when 
F = 0 or 1. In the limit F + —>• 0.5 then both observables 
vanish. It is interesting to note that a similar relation¬ 
ship between the two classes of D decays was found in 
Ref. [|] when considering the determination of the uni- 
tarity triangle angle 7 using —»• DK ± decays. 

Expressions (fl9l) and (1201) may be modified to allow for 
the possible contribution of direct CP violation. In this 
case the relations in Eq. (fTOl) no longer apply. Direct CP 
violation adds an additional magnitude and weak phase 
difference when considering the relations between the am¬ 
plitude of the D° and D a decay, and this additional mag¬ 
nitude and phase varies as a function of position in phase 
space. 


A 


eff 


(2Fl — l)y cos (f>cp (tcp - 

V t c 

(2 F' + - l)xsin 4>cp (tcp 4-) + 

V TCP J 


y A sin (j> C p r C p -+ 

TCP , 


zA cos <f>CP ( TCP + 


— ) 

TCP J 


( 21 ) 


where tcp and (f>cp are unchanged in their meaning 
and relate only to indirect CP violation, (2 F' + — 1) = 

eq \J K[ K 1 + c_i \J K- 1 K -1 and A = si y/Kjcl + 

s_i \JK- \ K- \. Hence the effect of the additional am¬ 
plitudes due to direct CP violation is contained within 
the terms F' + and A. In the limit of no direct CP viola¬ 
tion A —>• 0, and F' + —> F + . Since A must be small the 
third term in Eq. (1211) is negligible in comparison to the 
others. 

The expression for y^p becomes 


V e cp 


(2 F' + - l)y cos 4>cp (tcp + ~ — 
V TCP 

(2 F' + - l)x sin (f) C p (tcp 


TCP J 


yA sin 4> C p ( tcp H- — ) + 

TCP J 


X A COS 4>CP ( TCP — 


— ) 

TCP J 


( 22 ) 


Discussion and conclusions 

Measurements of Af? and y^ P performed with any 
self-conjugate multibody decay can be used to determine 
Ar and yep, respectively, provided that the CP content 
of the decay is known. The mode D n + n~n° is a 
very promising candidate for this purpose since the dilu¬ 
tion effects arising from the factor (2F + — 1) in Eqs. (fl9l) 
and © are < 10%, and it possesses a branching fraction 
that is around 3.5 times higher than that of D — » K + K~, 
the most common CP-eigenstate mode used for these 
measurements. Therefore this channel offers an oppor¬ 
tunity to improve the knowledge of Ar and yep signif¬ 
icantly, particularly at e + e _ experiments such as Belle- 
11, where the 7r° reconstruction efficiency is good. The 
relatively abundant four-body decay D —► 7r + 7r~7r + 7r - , 
which has a CP-even fraction of 0.737 ± 0.028 Q, also 
has the potential to be a high impact channel. The 
sensitivities of these channels are compared to those of 
the established CP-eigenstate decays, D —>• K + K~ and 
















4 


D —> 7 r + 7 r , in Table [T] assuming the same trigger and 
reconstruction efficiency for all. 


TABLE I: Relative uncertainties on Ar and yep for the multi¬ 
body modes D — > n + Tv~n° and D — > n + n~n + n^ compared 
with those of the CP eigenstate modes D —¥ K + K~ and 
D —>• 7r + 7r _ , assuming the measured central values of the 
branching fractions (BF) 0 and CP-even fractions (F+) 0. 
The uncertainties are all normalised to that of D —> K + K~. 



K+K~ 

7r + 7T 

+ 

1 

C 

1 

+■ 

1 

f 

TTF^xIcwT 

0.396 

0.1402 

1.43 

0.742 

F+ 

1 

1 

0.973 

0.737 

Uncertainty 

1 

1.68 

0.56 

1.54 


Another mode of potential interest is D —> Kg7r + 7r~7r°, 
which has a branching fraction of over 5% and comprises 
the CP-odd eigenstates Kg?] and Kgco as submodes, al¬ 
though its sensitivity cannot be assessed until its CP 
content is measured. This channel also has the feature 
of being Cabibbo favoured, which means that it is ex¬ 
tremely robust against any pollution from direct CP vi¬ 
olation. The extensively studied decay D —> /\g 7 r + 7 r~ 
is not suitable for an inclusive treatment, since it has a 
CP content of F + ~ 0.5, as is evident from examining 
the relative proportion of CP-even and CP-odd double- 
tagged events reported in a CLEO analysis performed to 
measure the Ci and s ? ; parameters [1_T|. 

The Belle collaboration has reported a model- 
dependent analysis of the mode D —> KgK + K~ that 
measures yep through comparing the CP-odd and CP- 
even regions of the Dalitz plot jl8(. Studies also exist 
that fit time-dependent amplitude models to the Dalitz 
plots of the decays D —> Kg 7r + 7r~ and D — > KgK + K~ in 
order to determine the mixing and CP-violation param¬ 
eters H-21|. Furthermore, proposals have been made 
of how to perform model-independent analyses of self¬ 
conjugate decays binned in phase space [22, 10. The 
method advocated in this paper is novel because it is in¬ 
clusive, model independent and suitable for those decays 
which are dominated by a single CP eigenstate, such 
as D — > . Inclusive analyses are experimentally 

more straightforward since there is no need to account 
for the position in phase space of each decay, provided 
that the acceptance is relatively uniform. 

As explained in Ref. fsj, self-conjugate multibody 
modes can also be used to measure the unitarity trian¬ 
gle angle 7 with B ± —> DK ± decays as long as F + is 
known for the mode under consideration. In cases where 
no measurement of F + exists from the charm threshold it 
is possible to obtain this information from a comparison 
of a measurement of y'fjp and the value of yep obtained 
from CP eigenstates, or indeed that of y itself, assuming 
negligible CP violation in the charm system. This strat¬ 
egy of using charm-mixing observables to help provide 
input for the 7 determination is similar to that already 


proposed for quasi-flavour specific states 


24]. 


In summary, inclusive measurements of the time evo¬ 
lution of mutibody self-conjugate charm decays offer the 
possibility to obtain significantly improved sensitivity to 
CP violation and mixing in the D°D° system. The ob¬ 
servables Ap ff and y^ p are simply related to those of the 
CP eigenstate case, Ar and yep, by a dilution factor 
(2P+ — 1), where F+ is the fractional CP-even content 
of the decay. This parameter may be measured in co¬ 
herently produced DD decays at the t/>(3770). One of 
the modes for which F + is known, D —» 7 r + 7 r _ 7 r°, has 
the potential to yield a more precise determination of 
Ap and yep than is possible with CP eigenstate decays, 
and another, D — > 7r + 7r _ 7r + 7r _ , also offers good sensi¬ 
tivity. Other promising channels exist with relatively 
high branching fractions and should also be exploited, 
provided that analyses at the ^(3770) show them to be 
dominated by a single CP eigenstate. Alternatively, mea¬ 
surements of y^p using these latter channels will allow 
their CP content to be determined, which is valuable in¬ 
put for the programme to measure the unitarity angle 7 . 
First results using this class of decays are eagerly awaited. 
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